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Abstract. Wc define parafcrmionic obscrvablcs in various lattice loop models, 
including examples where no Kramers- Wannier duality holds. For a particular 
rhombic embedding of the lattice in the plane and a value of the parafcrmionic 
spin these variables are discretely holomorphic (they satisfy a lattice version of 
the Cauchy-Riemann equations) as long as the Boltzmann weights satisfy certain 
linear constraints. In the cases considered, the weights then also satisfy the critical 
Yang-Baxter equations, with the spectral parameter being related linearly to the 
angle of the elementary rhombus. 



PACS numbers: 02.10.Ox, 05.50.+q, 11.25.Hf 

1. Introduction 

Holomorphic (and antiholomorphic) fields are the basic building blocks of conformal 
field theories (CFTs). They have simple short-distance expansions among themselves 
and with other fields of the theory. Although the prototype of such holomorphic field 
is the stress tensor T(z), which is present in all CFTs, many interesting CFTs contain 
holomorphic fields with fractional conformal spin. These are often referred to as 
parafcrmionic. The correlation functions of such holomorphic fields are necessarily 
power-behaved, and therefore, in local theories, can exist only in the massless, 
conformal case. 

Many CFTs are also believed to describe the scaling limit of critical lattice models. 
In order to understand the emergence of holomorphic fields in this limit, the simplest 
possibility is that there should exist analogues already in the discrete setting. These wc 
refer to as discretely holomorphic observables of the lattice model. A precise definition 
is given below. 

Recently a number of examples of such observables have been discovered in well- 
known lattice models. In models which possess a Kramers- Wannier duality symmetry, 
parafcrmionic observables can be defined in terms of suitable products of neighboring 
order and disorder variables [Tj. In the case of the Z^v, or clock, models, it was 
shown in [5] that such parafermions are indeed discretely holomorphic, but only at 
the integrable critical points identified by Fatcev and Zamolodchikov [3]. In this case 
the parafermions were defined algebraically directly in terms of the spin variables and 
Boltzmann weights of the model. 

A second example is the Q-state Potts model, which also possesses a duality 
symmetry. This model, in its Fortuin-Kastcleyn (FK) random cluster form, can also 
be represented in terms of a set of dense random curves [4] . In [5] , it was shown that 
the parafermions defined in terms of order-disorder products have support on these 
curves, and indeed may be considered as observables depending only on the curves. 
For example, the two-point function vanishes unless the two points happen to lie on the 
same curve. In [SJ[B] a simple argument was given that these parafcrmionic observables 
are discretely holomorphic, but only at the critical point of the Potts model (where it 
is also known to be integrable). 

In this paper we amplify these comments and enlarge the list of lattice models for 
which discretely holomorphic observables can be identified. In particular, wc consider 
various versions of the O(n) model which do not possess Kramers- Wannier duality and 
for which the definition of a parafcrmion in terms of an order-disorder product fails. 
However, these models all have a representation in terms of curves on the lattice, and 
using this we show that one may still define parafermionic observables, generalising 
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the construction in the Potts model. Moreover, in all cases, these turn out to be 
discretely holomorphic precisely when the model is integrable and critical. 

More generally, we may consider versions of these models where the weights are 
anisotropic. In this case, one would expect the scaling limit to be a rotationally 
invariant CFT only if the lattice is embedded in the plane in a particular way. 
In general, there is a one-parameter family of such incquivalent embeddings of a 
homogeneous lattice. For a given embedding, we find that discrete holomorphicity 
holds only for a particular set of weights. These weights, in all cases, satisfy the Yang- 
Baxter equation, with the usual spectral parameter being simply related to the angle 
of shear in the embedding. 

This is a remarkable result, since the requirement of discrete holomorphicity, 
as we will show, gives a set of linear equations in the weights for a fixed value of 
the spectral parameter. On the other hand, the Yang-Baxter equations are cubic 
functional relations between the weights at different values of the spectral parameter. 

Let us define more carefully the notion of discrete holomorphicity. Let Q be a 
planar graph (which will usually be a regular lattice) embedded in the complex plane, 
with vertices at points {zj}. Let F{z) be a function defined on the vertices of the 
medial lattice: the midpoints \{z% + Zj) of the edges (ij) of Q. Then F is discretely 
holomorphic on Q if 



where the sum is over the edges of each face, or plaquette, V of Q. Equation ([I]) 
may be thought of as a discrete version of Cauchy's theorem. However, note that, 
since there are fewer equations than unknowns, even in the discrete setting (flj is not 
sufficient to solve a boundary value problem on Q. (This should be compared with the 
example of the discrete Laplace equation.) In addition, if we try to take the scaling 
limit by covering the interior of some domain T> in the plane by a suitable sequence 
of graphs Q whose mesh size approaches zero, we may deduce that contour integrals 
vanish, but Morera's theorem implies that the scaling limit of F is analytic only if we 
first assume, or prove, that it is continuous. (In the Ising case, where the scaling limit 
has been proven [7], it was in fact necessary to adopt another less direct approach.) 

In the examples we discuss, the discrete holomorphicity follows entirely from 
local properties of the lattice model, and therefore we may in fact take F(z) to be 
any correlation function of the local parafermionic observable with other local fields at 
other locations, which may be thought of as a conditional expectation value. In what 
follows we shall not distinguish explicitly between the observable and its correlation 
functions. 

We note that Smirnov [6] has argued that if one can find suitable discretely 
holomorphic parafermionic observablcs of curves in lattice models, and also prove 
that they have a holomorphic scaling limit satisfying suitable boundary conditions, 
such that their dependence on z in simple domains is computable, then it follows that 
the scaling limit of the curves is Schramm-Loewner evolution (SLE K ) with a value of 
k related to the conformal spin. This programme has been carried to completion for 
the Ising model, for the curves forming boundaries of both the spin clusters and the 
FK clusters [7]. 

This application of discrete holomorphicity is not the purpose of the present 
paper. Rather it is to point out that it holds in a wider class of models than was 
observed up to now, and that it appears to be intimately related to integrability. 




(1) 
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However, we should note a general feature of the relation of our results to SLE. 
As was first observed in [5J [!|] , the statement that a curve starting on the boundary 
of a domain is SLE is equivalent in CFT to the assertion that the boundary operator 
which inserts the curve corresponds to a Virasoro highest weight representation with 
a null state at level 2, usually denoted by 02. l (or (pi, 2)- It has conformal weight 
^2,1 = (6 — k)/(2k). Now one of the important properties of a holomorphic bulk 
conformal field ip s (z) with spin s is that, when taken to the boundary, it gives a 
boundary field with the same total conformal weight, that is s. If the bulk field is an 
observable which depends on N curves, which meet locally at the point z but begin at 
different points {xj} on the boundary, this implies that the CFT correlation function 

N 

3=1 

is non-zero. By the fusion rules of boundary CFT, this implies that ip s must be a 
4>n +1,1 field, and therefore that 

N(2N + 4-n) 
s = /ijv+1.1 — o • 

ZK 

Thus if one knows the value of s one may infer the value of k. In this paper we give 
examples with N = 1,2, as well as cases when the curve cannot be simple SLE. 

We note that in this paper we shall consider homogeneous, translationally 
invariant lattices, but many of the arguments extend to so-called Baxter lattices, 
which may be embedded in the plane isoradially (that is all their faces are rhombi). 
These have been discussed with respect to the Zjv models in [2] and the Potts model 
in [5] . Discretely holomorphic functions in the Ising model are also studied in [TU] . 

The layout of this paper is as follows. In Section[2]we review the arguments for the 
FK representation of the Potts model. In Section[3]we give a parafermionic observable 
for Nienhuis' 0(n) model on the square lattice [TT], and show that it is discretely 
holomorphic precisely on the integrable manifolds [TT] . In Section 2] we extend this 
to a model Q~2] with two different kinds of loops which may cross each other. Our 
arguments in these last two cases are generalisations of those of Smirnov |13j for the 
O(n) model on the honeycomb lattice. 

2. The self-dual Potts model on the square lattice 

The Q-state Potts model is a classical spin model on the lattice. Each vertex carries 
a spin variable Si <E {1,2,..., Q}, and the Boltzmann weight for a spin configuration 
is 

W[{S i }]=l[[exp[J ij SiS^Sj)] , (2) 

(id) 

where the product is on the edges of the lattice, and Jij is the coupling constant 
attached to edge We consider the case of the square lattice C, with coupling 

constants Ji, J2 on the horizontal and vertical edges, respectively. The model can be 
reformulated [1] as a dense loop model on the medial lattice A4: this is the lattice 
consisting of the midpoints of the original lattice. The elementary configurations of 
the loop model are defined on square plaquettes V, which are the faces of the covering 
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Figure 1. Left: The original square lattice C of the Potts spins (dotted lines) 
and the covering lattice M* (full lines). Right: the deformed lattices, defined by 
the angle a. 





Figure 2. The elementary plaquette configurations for the loop model, and their 
Boltzmann weights. Any closed loop has weight \fQ. 



lattice Ai* , the dual of M. (sec left-hand side of Figure [I}. The Boltzmann weights of 
the elementary plaquettes (see Figure [5]) are (a r ,6 r ), where 

bi _ e Jl - 1 h VQ 

ax VQ ' a 2 e-*»-l' lj 
and r = 1, 2 respectively if the plaquette contains a horizontal or a vertical edge of the 
original lattice. Each closed loop carries a weight \J~Q. Equivalcntly, the loop model 
Boltzmann weights can be encoded in the i?-matrix R(u) = a(u) 1 + b(u) E, where E 
(the Temperley-Lieb generator) is the operator acting on loop diagrams as depicted 
on the right of Figure [2 and u is called the spectral parameter. 
In this paper, we will restrict ourselves to the self-dual line: 

(e Jl -l)( e j2 -l) = Q, Ji,J 2 >0, (4) 

in the case that the loop model is homogeneous. For < Q < 4, the Potts model has 
a second-order transition on this line. In this regime, we can write: 

v/Q = 2cos 7 , 0< 7 <|. (5) 

In [3 [6], a lattice holomorphic observable F s (z) was identified in this loop 
formulation. The observable F s (z) is defined on the midpoints of the plaquette edges 
as follows: 

F s (z)= J2 P(G)e- ise ^, (6) 
Ger(o.z) 
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where P(G) is the probability of the graph G, T(0, z) is the set of loop configurations 
for which the points and z belong to the same loop, and 8(z) is the winding angle 
of this loop from to z. More precisely, we fix arbitrarily a point on the edge 
of a plaqucttc (this could be in the bulk or on the boundary of the domain), and a 
direction for the loop segment passing through 0. The winding angle 8(z) is defined 
incrementally, setting 8(0) = 0, and adding ±^ for each elementary plaquette on 
the loop, according to the turn made by the loop. The real number s is the spin of 
the parafermion F s (z). In [5], the holomorphicity relation |T]) was established for F s 
defined in ©, at the isotropic self-dual point: J\ = J 2 = ln(l + y/Q). This is only 
true if the spin s satisfies: 

v/Q = 2sin^. (7) 

In this section, we adapt these results to the whole (anisotropic) self-dual line. 

On the one hand, since the loop plaqucttes obey a Tcmpcrlcy-Licb algebra, the 
i?-matrix satisfies the Yang-Baxter equations: 

Ri 2 (u)R 23 (v - u)R 12 (v) = R 23 (v)R 12 (v - u)R 23 (u) (8) 

for the parameterisation a(u) = sinu, b(u) = sin(7 — u). On the other hand, let 
us consider the following deformed model : this is the self-dual Potts model with 
(@| and anisotropic weights J\ ^ J 2 , defined on the rectangular lattice with vertical 
edges scaled by a factor cotan^ (see Figured]). Various arguments (see e.g. pHl UH\ ) 
support the statement that, in the continuum limit, the deformed model has the same 
behaviour as the isotropic self-dual Potts model, if the angle a is related correctly to 
the ratio 3 2 j J\. We shall obtain the following results for the deformed model: 

(i) F s is holomorphic on the lattice if the spin s satisfies the relation ||7J). 

(ii) There exists a linear relation between the spectral parameter u and the angle a. 

Let us describe in more details the deformed geometry: the plaquettes become 
rhombi with internal angles a, ir — a, and the increments of 8(r) are now ±a, ±(tt — a) 
on each plaquette. The left-hand side of ((T|) has only contributions from the loop 
configurations in which a closed loop C connects to two edges of the plaquette V. 
Let G be one such configuration, where the plaquette V is in configuration (a) of 
Figure [21 and G' the configuration differing from G only by the configuration of the 
plaqucttc V . Our method is to show that the contributions to the left-hand side 
of |T]) from G, G' cancel each other. There are essentially two inequivalent external 
connectivities of the loop outside of V 1 as shown in Figure [3] In the case (1), the 
probabilities of G, G' are related by P(G')/b = P(G)/(a^/Q), whereas in case (2), one 
has P(G')/ (by/Q) = P(G)/a. The requirement of cancellation of contributions from 
G, G' yields a linear system for the Boltzmann weights a, b: 

(1 + H)VQ a + (1 + n - A - fiX- 1 ) 6 = 

(9) 

(1+ (i— A" 1 - ^A- 1 ) a + (1 - fi\ )y/Q b = 0, 

where we have set A = e ms , fi = e 1Q ' s+1 '. The determinant of the system reads: 

A" x (l + - A^- 1 ) [A 2 + (Q - 2)A + 1] . (10) 

This has to vanish for the system to have a non-trivial solution. The last factor 
vanishes for A = — e ±2 ' 7 . We choose the minus sign, corresponding to the exponent: 

, 2 7 
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Note that this value satisfies the condition (J7J). The relative Boltzmann weight is then: 

b cos^ , 

(12) 



a cos 



Comparing with the parameterisation a{u), b(u) for the i?-matrix, we get: 
(s + l)a 7r 



(13) 



2 2 

Therefore, we recover the result of [15], that the angle a is linearly related to the 
spectral parameter u. Note that Q13p is different from the usual relation a = ttu/j. 
This can be corrected by setting the loop's angle increments on a plaquette to j3, tt — /3. 
Then the system ([5]) still holds, if we set (j, = e i ( Q+/3s ). The angle (3 can then be tuned 
to recover the appropriate relation between a and u. 

The scaling limit of the self-dual Potts model is described by a Coulomb gas with 
central charge and conformal weights given by: 

6(1 -9? 



c= 1- 



9 

(gr r'f (1 - gf 



(14) 



4.9 

where the coupling constant g is given by g = 1 — -. One can then write (jTTJ) as 
s = /i3,i- According to our earlier observation, this is consistent with the fact that 
there are two curves meeting at the observation point z. 



3. The O(n) loop model on the square lattice 

In this section, we find solutions of the holomorphicity equations ([1]) for the O(n) loop 
model on the square lattice [TT]. This is a dilute loop model, defined by the vertices 
in Figure [4] where each closed loop carries a weight given by: 

n = — 2 cos 2rj , < r\ < — . 

Note that we have grouped the anisotropic weights so that there is symmetry under 
reflections in the diagonal axes. This symmetry is preserved when the plaquettes 
are deformed into rhombi. Since every loop configuration has an even number of 
plaquettes of type U\ or u 2 , the change (1*1,^2) ~~ * (~ u i,~ U2) does not affect the 
Boltzmann weights. 

We consider the observable F s (z), defined similarly to ©, except that the sum 
is now over graphs where there is an oriented open path going from to z. A similar 
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Figure 4. Vertices of the O(n) loop model on the square lattice. Each closed 
loop has a weight n. 




t T T T 

Figure 5. Loop configurations with one edge of the plaquette V connected to 
point 0. 



observable has been considered by Smirnov [T3] for the case of the hexagonal lattice. 
Now consider the contributions to the holomorphicity equation ([T]) where the first time 
the oriented curve, starting at 0, enters the chosen plaquette V is through a particular 
edge, for example the lowermost in Figure 0- There are four inequivalcnt external 
connectivities as shown. We imagine these to be fixed, summing over all internal 
configurations consistent with them. This yields the linear system for the Boltzmann 
weights: 

t + fiui — fi\^ 1 u2 — v = (15a) 

— A _1 ui + i%U2 + X/iv — fj,X (w-y + nw?) = (15&) 

nu\ — Xii2 — yuA~ 2 v + n(nw\ + W2) = (15c) 

— /iA~ 2 «i + /iA«2 + nv — A~ 2 wi — A 2 W2 = , (15rf) 

where we have set: A = e 1TS , ip = (s + l)a , /i = e iv . For real Boltzmann weights, (|15al - 
Il5a| are four complex linear equations for six real unknowns (t, ui, u%, v, Wi, W2), and 
we have the relations: 

Im [(n + 1) CEil - X^ 1 Cig) + n- 1 (HHcJ)] = 

Im [Apr 1 (A 2 - nX- 2 ) (HHJ) + ^{nX 2 - A" 2 ) dUcJ) - {n 2 - 1) (ITBUl ] = . 

Thus, we can generally reduce (I15alll5dl) to a 6 x 6 real system. 

There are two classes of solutions, for vanishing and non-vanishing v. First, if 
v = 0, then the configurations corresponding to (Il5c£j) never occur, and so this equation 
does not hold. In the special case n = 1, there exists a non-trival solution for any 
value of s: 

t — sin7rs , u\ = sin(cp — 7rs) , 1*2 = simp , Wi + u>2 = sin7rs . (16) 

This model can be mapped onto the six- vertex model (see Figure 15]), with weights 
u>i = L02 = sin(p — tts),u>3 = L04 = sine/?, W5 = u>e = sin7rs. The corresponding 
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Figure 6. Mapping of the O(n) model onto the six-vertex model for n = 1, V = 0. 



anisotropy parameter [4] is A = cos7rs. This is an example of a model admitting 
a holomorphic observable on the lattice, but for which the scaling limit of the 
corresponding curve cannot be described by simple SLE. This is because, for ordinary 
SLE, the central charge of the CFT is directly related to the SLE parameter k [HI E] 
and hence to the conformal spin s: c = 2s(5 — 8s) /(2s + 1). In the present case, since 
the boundary conditions for the six-vertex model are not twisted, its scaling limit has 
central charge c = 1 for all A. However the conformal spin s varies continuously with 
A. Therefore the scaling limit of the curve can be SLE, with k = 4, for at most 
one value (in fact A = l/\/2-) We conjecture that other values of A correspond to 
SLE(4,p). 

For v = 0,n 7^ —1, we get a 5 x 5 linear system, with determinant (n 2 — 
l) 2 smtp sin((p — its). Imposing sinip = yields (s + l)a = mir and in turn s = to', 
where to, to' are integers. Thus, for the solution to exist at any value of a, we have 
to set s = — 1. The Boltzmann weights are then: 

t = -ui - u 2 , w\ = —u\ , w 2 = —U2 ■ (17) 

The solution of the case s'm(ip — its) = is similar, and leads to the same Boltzmann 
weights and spin s = — 1. If we change the sign of Ui,u 2 , then the model (|17[) is 
equivalent to the dense loop model of Section [51 with parameters ^/Q = n + l,a = 
Ui,b = u 2 . To see this, fill empty spaces with loops of weight 1 (ghost loops). 
The local weights do not depend on the type of loops involved (actual or ghost 
loops), so each loop has an overall weight n + 1. As a consequence, the dense loop 
model has a lattice antiholomorphic observable (s < 0), besides the holomorphic 
one found in [5J. However, several arguments rule out the hypothesis that F s= ^i(z) 
corresponds to an antiholomorphic field in the continuum limit. First, F s= ^i(z) 
is lattice antiholomorphic for any Q > 0, whereas it is well known that the self- 
dual Potts model is only critical for < Q < 4. Furthermore, the ratio u\ju 2 
in (|17p does not depend on the angle a, which means that the same model has 
an antiholomorphic observable for any deformation angle: this is not acceptable 
physically in the continuum limit. So we conclude that, in the case v = and generic 
n 7^ —1, the holomorphicity conditions ([T]) for the dilute O(n) model merely lead to 
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the case of the dense loop model, but the corresponding F s {z) is not a candidate for 
an antiholomorphic field in the continuum limit. 

Let us now discuss the solutions of second class (w ^ 0), for a generic value of n. 
We get the 6x6 real system: 

{Re (fTgqj) . Re (fT56)) . Im (TT5b|) , Re (fT^j) . Im (|T5^|) . Re (fT5d|) } , 

with determinant: (n 2 — 1) sinos sin(o3 — 7rs) (2 cos47rs — 3n + ?i 3 ) . Non-trival solutions 
exist if the spin satisfies: 

cos Airs = cos 677 . (18) 

The various solutions to (fT5)) can be parameterised by extending the range of 77 to 
[— 7T, 7r], and setting: 

Then, we get the second class of solutions, with Boltzmann weights: 

t = — sin(2o3 — 377/2) + sin5?7/2 — sin3?//2 + sinry/2 

u\ = —2 sin 77 cos (3r)/2 — tp) 

Ui = —2 sin 77 sin 03 

7; = —2 sin 03 cos (077/2 — tp) 

W\ = —2 sin(o3 — ?/) cos (3t//2 — tp) 

W2 — 2 cos (77/2 — 03) sin 03 . 

A remarkable fact is that the weights ([20]) are a solution of the Yang-Baxter equations 
for the O(n) loop model on the square lattice. Indeed, after a change of variables 
03 — > i/'+( 7r + ? ?)/4, (|2"0"1) coincides with the integrable weights in [TT]. So, by solving the 
holomorphicity equations (|15alll5djl on a deformed lattice, we recovered the integrable 
weights. 

Now we interpret our findings in terms of CFT. First, in the regime < 77 < ir 
(branches 1 and 2 in |16j ) , on the basis of numerical diagonalisation of the transfer- 
matrix [T^] and analysis of the Bethe Ansatz equations [TTlll8| . it has been argued that 
the model (|2"0|) has the same continuum limit as the 0(n) model on the hexagonal 
lattice |19| : it is described by a Coulomb gas of coupling constant g = 277/71", with 
central charge and critical exponents given by (|14|) . The value of the spin can then 
be written as s = /i2,i- In the regime — tt < 77 < 0, the continuum limit is no longer 
described by a simple Coulomb gas. Rather, it has central charge: 

3 6(1 -g'f 

C= 2~^^> (21) 

where g' = 2(ir + ij)/tt. The relation between s (fT9|) and c (f2Tj) . which is c = 
(2 — 5s — 16s 2 )/(2s + 2), is not consistent with SLE, so this is another example of a 
model with a holomorphic observable whose curves are not described by simple SLE. 
Note that in this case c may be greater than one, so some variant of SLE such as that 
considered by Bettelheim et al. [20] may be the correct description. 
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Figure 7. Vertices of the loop model on the square lattice. Each closed 

loop has a weight n. 



4. The C { 2 ] loop model 

A simple generalisation of the loop model of Section[5]to a model with two loop colours 
was introduced in [12] . It is a dense loop model on the square lattice, where each loop 
can be either black or grey, and carries a weight n given by: 

n = — 2 cos 2i] , < 77 < — . 

The vertices of the model are shown in Figure [7J For fixed n, the model has a one- 
parameter family of Boltzmann weights that satisfy the Yang-Baxter equations [T2"] . 

To describe the observable F s (z) that will be shown to be holomorphic, we first 
introduce a special kind of defect in the loop model: this defect sits on the midpoint 
of a plaquette edge, and consists in a change of colour for the loop passing through. 
The observable F s (z) is then defined similarly to ([6]), except that the sum is now over 
graphs G where and z belong to the same loop, and each of these points carries a 
defect. The angle 6(z) is defined as the sum of the windings of each loop strand from 
to z. For a given plaquette V in fl}, the graphs that contribute are those for which 
two edges of the plaquette V are connected to the point 0, and one edge contains a 
defect. The holomorphicity conditions (fTJ) yield the system: 

nu\ — \ii2 — fi\~ 1 v + fi(nwi + u'2) = (22a) 
— A _1 ui + nu2 + \iv — /jA _1 (wi + nw^) = (226) 
— /iA _1 ui + fxu% + nv — A _1 wi — Au>2 = , (22c) 

where we have set: A = e 2l1TS , ip = (2s + l)a , [i = e llf . For real Boltzmann weights, 
we have the linear relation: 

Im [M -1 (nA - A" 1 ) (J22it) + A/J _1 (A - nA" 1 ) ([226l + (n 2 - 1) = . 

Thus, for n 2 ^ 1, we can ignore Im (|22cp . and the remaining 5x5 real system has 
determinant: (n 2 — l)siny sin(y> — 2irs) (2cos47rs — in + ?i 3 ). It has a non-trivial 
solution for: 



cos Airs = cos 677 



(23) 
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There are two branches of solutions, which can be obtained by setting s = (3t? — 
7r)/ (27r), and extending the range of rj to [0, n]. The Boltzmann weights then read: 

ui = sinjysin^ — 3rj) 
«2 = — sin 77 sin (p 

v = — sin<psin(<^ — 377) (24) 

wi = — sm(ip — 77) siii(<p — 3??) 

W2 = — sin(<p — 277) sin tp . 
These are exactly the integrable weights given in [T^] , with the change of variables 
(77 — * A, if — > u). So the C2 1 00 P model is another example where the solution of the 
lattice holomorphicity equations fl} satisfies also the Yang-Baxter equations. 

5. Summary 

In this paper we have identified several more examples of discretely holomorphic 
parafermionic obscrvables in lattice models, and shown that this requirement always 
appears to pick out the critical points which are also integrable in the sense of Yang- 
Baxter. It would of course be important to have a more general understanding of this. 
It may give a simpler route to finding new integrable models. These parafermions 
are natural discrete candidates for parafermionic holomorphic conformal fields in the 
scaling limit, and thus can also suggest new structures in particular CFTs. In some 
cases the models we have considered correspond to CFTs with c > 1, so the scaling 
limit of the corresponding curves is probably described by some modification of simple 
SLE. The loop models we have considered may also be mapped to generalised restricted 
solid-on-solid (RSOS) models, and, from this point of view, it would be interesting to 
understand the coset construction of the corresponding CFT and its relation to the 
parafermionic fields. 
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